Epitaxial growth under oblique incidence 
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A continuum one dimensional model of homoepitaxial growth under oblique incidence is inves- 
tigated. We carried out numerical integration of a continuum equation incorporating a shadowing 
search algorithm. The interplay between the Ehrlich-Schwoebel (ES) effect and shadowing is clearly 
highlighted. It was found that different growth phases are separated by a well defined crossover 
time after which the ES and shadowing are the dominating mechanisms. Also, the model identified 
the existence of a transition period where deep grooves develop on the surface before the shadowing 
regime is fully developed. We found that growth under oblique angles accelerates the coarsening of 
mounds on the surface and we determined the corresponding critical exponent. 

PACS numbers: 



Epitaxial growth is a process where atoms are de- 
posited on a crystaUine substrate to undergo diffusion 
on terraces or at the edge of atomic steps as well as at- 
tachment and detachment from step edges. These micro- 
scopic processes determine the surface morphology of the 
growing film[ll,|3. It is well known that the formation of 
mounds during Homoepitaxial growth is a result of asym- 
metric attachment-detachment kinetics i.e. at the edge 
of an atomic step, atoms have to overcome an energy 
barrier to move downward. This energy barrier is called 
the Ehrlich-Schwoebel (ES) barrierj^ In spite of the 
fact that many experimental observations of the surface 
morphology of films grown by molecular beam epitaxy 
have been explained by the above relaxation processes, 
it was recently reported that geometrical effects can play 
an important role when growth is performed under 

oblique incidence. In Cu/Cu(100) epitaxy for example, 
it was reported that large incidence angles 9 of incom- 
ing atomic flux result in the formation of asymmetric 
mounds for 9 = 70° and asymmetric ripples on the sur- 
face for 9 = 80° [5:]. In addition, the crystallographic 
orientation of the mounds / ripples facets depend strongly 
on temperature. These observations were confirmed by 
computer simulations which reveal that they are purely 
determined by geometrical effects 6]. 

Many theoretical aspects of growth in which shadowing 
plays an important role were introduced through contin- 
uum models0,@]. An example, is the sputter deposition 
of films where the observed surface morphologies were 
attributed to the influence of shadowing which induces 
an instability in the growth process fs*]. A simple contin- 
uum model that describes this effect is proposed in Q, 
in which the smoothening effect of diffusion is included 
as well. In this model, the height h of the surface grows 
according to the following equation: 



dh 
'dt 



^Ky^h + F0 + r]{x,t) 



(1) 



process respectively, while 77 is Gaussian white noise sat- 
isfying < ri{x,t)ri{x',t') >= D5{x - x')5[t - t'). The 
angle 9[x, {h}) is the exposure angle which depends on 
the position and the entire surface height h. In this exam- 
ple, the resulting surface has a columnar structure with 
deep grooves between neighboring columns, with a well 
defined characteristic thickness. The tops of the columns 
are growing with a constant speed along the surface nor- 
mal. 

Despite these efforts, many aspects of epitaxial growth 
under oblique incidence remain unclear. Furthermore, 
the shadowing effects have been largely unexplored ex- 
perimentally due to a lack of theoretical predictions of 
the experimentally observable quantities that discrimi- 
nate different growth regimes. Here we propose a one di- 
mensional continuum model of homoepitaxial growth un- 
der oblique incidence that incorporates the ES effect, dif- 
fusion and shadowing. We clearly uncover the interplay 
between these effects. We show a clear distinction be- 
tween the linear regime and the non-linear regime where 
the shadowing becomes significant, through the identifi- 
cation of a well defined crossover time between the two 
regimes. We identify the existence of a transition period 
where deep grooves develop on the surface before the 
shadowing regime is fully developed. Finally, we found 
that growth at oblique incidence accelerates the coars- 
ening of mounds on the surface and we thereafter deter- 
mined the corresponding critical exponent. 
A phenomenological continuum model describing the sur- 
face growth incorporating the ES effect, diffusion and 
shadowing can be formulated in one dimension by: 



dh 
'dt 



-7^ = - Vjd - Vj, + 3(9, X, h, m)F + r^{x, t) (2) 



Model ([T]) assumes low vapor pressure and ignores local 
non-linearity caused by lateral growth. The parameters 
F and K represent the deposition rate and the diffusion 



where h is the single valued surface height, jd is the ES 
destabilizing current, js is Mullins stabilizing diffusion 
current, F is the deposition rate and m = dh/dx is the 
slope. A model for the currents ja and js can be ex- 
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pressed as0, [ol. Hoj: 
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Here, v and K are positive constants, and the current 
function / satisfies f{x) = 1 — . Without the shadow- 
ing term equation Q predicts a mound-hke profile where 
the symmetric mounds are formed during the initial stage 
of the growth as a result of the competition between the 
ES effect and surface diffusion. Later in time, coarsening 
of mounds occurs with slope selection i.e. the slopes of 
mounds converge to a single value TOq @, • 
The function S' is a non-local shadowing term which in 
general depends on the incidence angle 9 (which is the 
angle between the incoming flux and the y-axis) of the 
incoming flux, the surface height and the surface slope 
at a given position on the surface profile. It can be de- 
fined as the probability that a point M(x,h,m) of the 
evolving profile receives the incoming flux which hits the 
surface under oblique incidence. This function arises in 
the field of optics and was extensively studied to take 
into account the shadow effect during the scattering of 
electromagnetic waves on randomly rough surfaces ( see 
for example [12] and references therein). The general ex- 
pression for the shadowing function S is given bvjT3|: 



S{e, M) = iJ (cot 6* - to) exp 



f{e\M,x)dx\ (4) 



where H is the Heaviside function i.e. H{x) = for 
X < and H{x) — 1 for x > 0, and M is a point on the 
profile having the height h and the slope m. The function 
f{9\M, x) is the conditional probability that the incoming 
flux intersects the surface in the interval [x,x + dx] with 
the knowledge that it does not hit the surface in the 
interval [0,a;]. This function is expressed as: 



fiO\M,x) 



{m! ~ ^)p{h' ]m'\h]m)dm' (5) 



where = cot 9 and h' = h + j^ix. The function p is the 
joint probability distribution of the heights {h, h') and 
the slopes (to, to'). Analytical expressions of S are possi- 
ble to derive for static rough surfaces; for example Smith 
and Wagner 14 1 and recently Bourilier and Berginc [T§| 
derived S for Gaussian one dimensional profiles. 
For evolving profiles, analytical expressions for S are in- 
tractable if not impossible to obtain, especially when non- 
linearities are present. In the latter case the joint proba- 
bility in ([5]) may not be the product of individual height 
and slope probabilities. It is however possible to com- 
pute numerically and with great accuracy this function 
for any profile, using a numerical algorithm [l5]. This 
purely geometric algorithm determines if a point in the 
profile is shadowed or illuminated. We use this algorithm 
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FIG. 1: Top plot: evolution of the quantity v =< h > —Ft 
as a function of time for 9 = 80°, 70°, 60° and 55°, and for 
1/ = K = 0.2. Bottom plot: the crossover time t* as a function 
if the incident angle showing the logarithmic dependence of 
t* on fi — cotO. 



to determine S' in ^ . In our model, values of S' = 1 are 
attributed to illuminated areas while values of 5 = are 
attributed to shadowed areas. 

We integrated equation ^ numerically using Adams- 
Bashforth scheme [l6| and imposed periodic boundary 
conditions. At each time iteration, the function S is de- 
termined at each point of the profile, before the height 
is updated. The initial surface is perturbed with a white 
noise. The numerical integration of Q revealed the ex- 
istence of an incident angle 9* below which the shad- 
owing effect is irrelevant. This critical angle is mainly 
determined in the linear regime and is dependent on the 
ratio v/K. The smaller this ratio, the larger is the crit- 
ical angle. Typically, 9* is in the range 45°-50° for v/K 
varying from 10 to 1. To monitor the effect of oblique in- 
cidence on growth we followed the evolution of the mean 
height of the evolving profile, a quantity which is ac- 
cessible experimentally with the help of scanning probe 
microscopes. Figure [T] shows the plot of w =< h >- 
Ft versus time( here < /i > is the mean height), for 
9 = 80°, 70°, 60° and 55°, and for u = K ^ Q.2 and 
F ~ 1. For incident angles larger than the critical an- 
gle, one can distinguish two regimes: the first one where 
< h >= Ft corresponding to the linear regime; the sec- 
ond corresponding to a growth phase where the shad- 
owing influence becomes relevant. This clear distinction 



identifies a time t* separating the two regimes. For an- 
gles smaller than the critical angle, the time t* becomes 
extremely long and the separation of the linear and the 
shadowing regimes is no longer clear due to the domi- 
nance of the ES non-linearities. In this case, the growth 
of the surface profile is mainly determined by diffusion 
and the ES currents since equation ([2]) is reduced to the 
well known MBE equation. The time t* can be estimated 
from geometrical considerations as follows. During the 
linear regime the typical height is a which is the surface 
width defined as a = Si 

where I is the typical distance between mounds given 
by I = 2TTy/{2K/i^y2]. If we consider two neighboring 
mounds of heights hi-i and hi (i is the position on the 
profile) then the mound {i — 1) shadows the mound i 



when cot9 < 



hi — I— hi 



dh 
dx 



t > t* 



^ a/l. This gives: 
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—\n[cot9) 

V 



(6) 



Figure [T]-b shows the logarithmic dependence of t* on 
/i — cot6. The time t* was identified as the time when 
the quantity v =< h > —Ft drops from zero to negative 
values as clearly shown in figure [TJ 

Beyond the time t*, nonlinearities caused by shadow- 
ing and the ES barrier become significant. Straight after 
the linear regime and before non-linearities become fully 
developed, the height h goes through a phase where deep 
grooves form. This is because valleys become shadowed 
and stop growing whilst the mound tips continue to grow. 
At this stage, the ES effect is still weak. As soon as the 
ES effect starts to be relevant, the deep grooves acquire 
a mean velocity due to the ES non-linearities, and the 
grooves morphology disappear, resulting in a morphology 
dominated by asymmetric mounds or columns at glanc- 
ing angles. Figure [2] shows an example of the evolution of 
the profile's morphology, demonstrating the three phases 
for the following parameters : 9 = 80°, v = 0.1, K = 0.2 
and F = 1. This morphological behavior is true for any 
9 > 9*. 

These observations can be quantified by considering 
the surface width a. In figure [3] the time evolution of 
this quantity is showing the clear distinction between 
the above mentioned growth phases. The parameters 
used here are = .1, K = .2 and F = 1. The early 
time growth is well predicted by the linear theory i.e. 
a = exp(i/t/Z^). The deep groove phase (indicated by the 
shaded area in figure [3]) induces a sudden increase of the 
surface width; this phase is followed by a phase where the 
surface width evolves following a power law e.g. cr ~ 
In the absence of shadowing, the scenario predicted by 
equation ([2]) is as follows: after the linear phase where 
the surface profile undergoes an exponential growth and 
where mounds form, a coarsening phase develops. During 
this phase, mounds coarsen and the typical mound size 
A follows a power law increase in time i.e. X ^ , where 
(3 = l/sQ. In the presence of shadowing, this scenario 




FIG. 2; Portion of the evolving profile for three different 
times: t=10, 20 and 100 for an incident angle of 6 = 80°. 
At t=10, the growth is dominated by the linear instability. 
The intermediate phase (t=20) is characterized by the devel- 
opment of deep grooves on the surface as the shadowing starts 
to become relevant. Later on (t=100), the deep grooves dis- 
appear and the surface settles to a morphology dominated 
by the ES effect and shadowing, characterized by develop- 
ing mounds and coarsening. The parameters used here are 
f = 0.2, K = 0.2 and F = 1. The system size was = 500. 



y~yln(t) 



Mounds/columns 




FIG. 3: Time evolution of surface width a showing the clear 
distinction between three growth phases, for an incident angle 
of S = 75° and a system size A'^ = 500. The parameters used 
here are = .1, K = .2 and F = 1. The early time growth is 
well predicted by the linear theory i.e. a ~ exp{ut/l^). The 
deep groove phase (indicated by the shaded area) induces a 
sudden increase of the surface width, followed by a growth 
phase where it evolves following a power law e.g. a ^ F . 
The inset shows a log-log plot of a vs t for long integration 
times. The power law exponent 7 was determined by non- 
linear regression fit giving a value 7 = 0.87 ± 0.001. The 
result was averaged over 20 realizations. 
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FIG. 4: Log- log plot of the mound size A versus time for 
incident angles 9 = 80°, 75°, 70°, 65°, 60° and 55°. The 
power law exponent was determined by regression fit which 
gives the value of /3 = 0.81 ± 0.001. 



remains similar. Indeed the coarsening is still persistent 
but evolves more rapidly than in the absence of shadow- 
ing. To show this, we performed long time integration 
of equation ^ and computed the height autocorrelation 
function C{r) —< h{x + r)h{x) The mounds lateral 
size A is given by the zero crossing of C (r) . Figure [3] 
displays a log-log plot of the mounds lateral size versus 
time for different incident angles. The size of the mounds 
follows a power law behavior A ~ t'^ with the exponent 
/3 = 0.81 ± 0.001. Notice that the value of the exponent 
is the same for all angles (larger than the critical angle). 
This tells us the coarsening is faster than in growth in 
normal incidence where the value of the exponent is 1/3, 
which is re-confirmed by the integration of equation Q 
for angles smaller than 6* . 

The draw back of the model ^ is it does not in- 
clude faceting to account for crystallographic effects and 
it is restricted to solid-on-solid (SOS) picture. Shim 
and Amar used kinetic Monte Carlo method which 
included faceting to simulate homoepitaxial growth of 
metal(lOO) under oblique incidence and found that many 
aspects of the surface morphology can be explained by 
purely geometrical effects induced by shadowing and 
demonstrated the appearance of ripples and rods with 



sides dominated by (lll)-oriented facets. Model ([2]), al- 
though it is one dimensional, qualitatively shares some 
features observed in these simulations; in particular the 
power law behavior of the surface width and the mounds 
size as well as the development of asymmetric mounds. 
In addition, model ^ provided an insight onto the in- 
terplay between different growth mechanisms involved in 
the growth process. Another advantage is that it can 
simulate long time behavior of the surface morphology 
without exorbitant computer resources. A natural pro- 
gression towards simulations which can directly be com- 
pared to experimental observations is to extend model 
([2]) to two dimensions. This will have an impact on the 
experimental design of naturally evolving nanostructures 
without resorting to expensive lithographic methods such 
as electron beam lithography. 
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